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Abstract 

High-sensitivity and liigli-througliput mutation detection techniques are useful for screening the homoplasmy or 
heteroplasmy status of mitochondrial DNA (mtDNA), but might be susceptible to interference from nuclear mitochondrial 
DNA sequences (NUIVlTs) co-amplified during polymerase chain reaction (PCR). In this study, we first evaluated the platform 
of SURVEYOR Nuclease digestion of heteroduplexed DNA followed by the detection of cleaved DNA by using the WAVE HS 
System (SN/WAVE-HS) for detecting human mtDNA variants and found that its performance was slightly better than that of 
denaturing high-performance liquid chromatography (DHPLC). The potential interference from co-amplified NUMTs on 
screening mtDNA heteroplasmy when using these 2 highly sensitive techniques was further examined by using 2 published 
primer sets containing a total of 65 primer pairs, which were originally designed to be used with one of the 2 techniques. 
We confirmed that 24 primer pairs could amplify NUIVlTs by conducting bioinformatic analysis and PCR with the DNA from 
143B-p° cells. Using mtDNA extracted from the mitochondria of human 143B cells and a cybrid line with the nuclear 
background of 143B-p° cells, we demonstrated that NUMTs could affect the patterns of chromatograms for cell DNA during 
SN-WAVE/HS analysis of mtDNA, leading to incorrect judgment of mtDNA homoplasmy or heteroplasmy status. However, 
we observed such interference only in 2 of 24 primer pairs selected, and did not observe such effects during DHPLC analysis. 
These results indicate that NUIVlTs can affect the screening of low-level mtDNA variants, but it might not be predicted by 
bioinformatic analysis or the amplification of DNA from 143B-p° cells. Therefore, using purified mtDNA from cultured cells 
with proven purity to evaluate the effects of NUMTs from a primer pair on mtDNA detection by using PCR-based high- 
sensitivity methods prior to the use of a primer pair in real studies would be a more practical strategy. 
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Introduction 

Human mitochondrial DNA (mtDNA) is a circular genome with 
16569 base pairs (bp), and exists as hundreds to thousands of 
copies in the mitochondrial matrix. It encodes genes for 2 rRNA 
genes, 22 tRNA genes, and 13 mRNA genes for subunits of 
Complexes I, III, IV, and V in oxidative phosphorylation. Several 
genetic features of mtDNA differ from those of nuclear DNA 
(nDNA), such as maternal inheritance, replicative segregation 
during cell division, and high sequence evolution rate [1]. 
Pathogenic mutations of mtDNA transmitted maternally can 
cause genetic diseases, whereas inherited ancient adaptive variants 
may increase the risk of certain diseases [2]. Somatic mtDNA 
mutations in aging-related diseases and cancers in humans have 
also been investigated extensively [2,3]. Because decision-making 
on homoplasmy or heteroplasmy status of mtDNA mutations can 
have implications, such as the determination of maternal 
inheritance or the differentiation between somatic mutations and 
genetic drift in cancer patients, detection of low-percentage 



heteroplasmy of mtDNA by high-sensitivity techniques, such as 
denaturing high-performance liquid chromatography (DHPLC) 
and pyrosequencing (PSQ), has become an increasing research 
interest, as we discussed previously [4]. PSQ^ is a powerful 
technique for the quantification of mtDNA heteroplasmy for a 
known single nucleotide polymorphism (SNP) or mutation without 
the need of establishing a standard curve [4,5], whereas allele- 
refractory mutation system (AMRS)-based quantitative polymer- 
ase chain reaction (qPCR) represents a highly sensitive technique 
that requires a standard curve made by mixing 2 DNA samples 
containing different alleles [6] . However, these kinds of methods 
are not suitable for the search on unknown variants or mutations, 
which are often present at low-level heteroplasmy as well, over 
entire mtDNA in various studies. Generating a standard curve for 
any potential novel SNP or mutation detected by direct 
sequencing under such scenario is also not practical. In contrast, 
DHPLC represents a kind of sensitive and high-throughput, 
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although quahtative, technique to screen heteroplasmic variants 
for entire mtDNA [7,8]. 

The detection of mutation by using the CEL 1 endonuclease 
from celery, which exerts 2 independent cuts at 3 '-side of 
mismatched nucleotides on 2 strands of heteroduplexed DNA 
without sequence specificity [9], has been commercialized as 
SURVEYOR Nuclease (SN) in die SURVEYOR Mutation 
Detection Kit. The platform of SN digestion followed by the 
detection of DNA fragments by using the WAVE HS System, 
designated as SN/WAVE-HS in this paper, is a high-sensitivity 
detection technique, in which an injector for the fluorescent dye 
and a fluorescent detector are equipped on the WAVE System 
used in traditional DHPLC analysis. Janne et al. has applied this 
platform for high-sensitivity mutation screening of the epidermal 
growth factor receptor (EGRF) gene in human cancer specimens 
[10], but it has not been applied for the detection of mtDNA 
mutation. Although Bannwarth et al. have developed 1 7 primer 
pairs for screening mutations of entire human mtDNA by using 
SN analysis, they detected cleaved DNA fragments by using 
agarose gel electrophoresis [11,12]. Because SN/WAVE-HS 
analysis does not require the testing on optimal temperatures 
and can be used to analyze amplicons much longer than 500 bp, it 
is very likely to be a more superior technique than DHPLC for the 
efficient screening of unknown heteroplasmic variants or muta- 
tions in entire mtDNA. 

The abbreviation NUMTs (or Numts) can refer to nuclear 
mtDNA sequences or mt DNA-like sequences in the nucleus [13], 
nuclear mitochondrial pseudogenes [14], or nuclear DNA 
sequences of mitochondrial origin [15]. Studies have indicated 
that NUMTs present in modern eukaryotes originated from the 
transfer and integration of ancient mtDNA fragments into nDNA 
at different stages during evolution and therefore have similar, but 
not identical, sequences to modern mtDNA sequences because of 
the much faster rate of mutation of mtDNA compared with nDNA 
[14]. According to results of bioinformatic analyses in different 
studies, human NUMTs, estimated to occupy 0.016% of nDNA 
[16], with the lengths ranging from 28 to 14 915 bp [13], exist in 
all human chromosomes as a single copy, multiple copies, or in 
highly rearranged forms, and match with all regions of human 
mtDNA sequences at various frequencies and with different 
similarities [13,17,18]. However, results from bioinformatic 
analysis of human NUMTs, such as lengths of NUMTs and 
location on chromosomes, varied among the diflerent studies 
because of the use of different searching methods or matching 
parameters. 

The possible misinterpretation of novel heteroplasmic mtDNA 
mutations in human diseases because of NUMTs [14,19,20] and 
the development of strategies to avoid interference from NUMTs 
during mtDNA detection [21,22] have recendy attracted consid- 
erable research attention. However, certain assumptions and 
approaches employed in previous studies require further exami- 
nation. First, the disputation of the novelty of a reported mtDNA 
variant based on the results of database searches, Internet 
searches, and phylogenetic analyses might not necessarily indicate 
that the variant was indeed from NUMTs unless otherwise proven 
by sohd experimental data. Because the copy number of the 
mtDNA outnumbers that of NUMTs on nDNA in regular DNA 
samples and direct sequencing has poor sensitivity for the 
detection of heteroplasmy, matching hc-tcroplasmic variants from 
sequencing results to predicted NUMTs by using bioinformatics 
approaches might not disprove the presence of actual mutations. 
Second, because the same amount of DNA from p" cells has much 
higher absolute amount of nDNA than the nDNA in DNA 
samples from regular cells containing both nDNA and mtDNA, 



the amplification of NUMTs from the DNA of p cells by a primer 
pair should not automatically infer the interference of NUMTs on 
mtDNA detection by any technique using this primer pair. Third, 
because some regions of mtDNA display high similarity to 
predicted NUMTs, the design of a usable primer pair specific to 
mtDNA might be problematic if low matching with NUMTs 
databases and inability to amplify the DNA of cells are the 
prerequisites. Finally, the use of pure mtDNA without nDNA 
contamination should be able to prove or disprove the possible 
interference of NUMTs on mtDNA detection by using a PCR- 
based technique. However, it is essential to demonstrate the purity 
of such mtDNA fractions through well-controlled PGR examina- 
tions because any nDNA contamination detectable by the PGR is 
likely to be detected by high-sensitivity techniques. 

In this study, we first established our SN/WAVE-HS methods 
for mtDNA detection, and compared the performance of SN/ 
WAVE-HS analysis with DHPLG analysis for the detection of 
different percentages of SNP by using 2 cloned plasmid DNA 
samples containing mtDNA fragments from human blood DNA. 
We then examined possible interference from NUMTs on the 
interpretation of heteroplasmy or homoplasmy during mtDNA 
detection by the DHPLG and SN/WAVE-HS techniques. Two 
published primer sets originally designed for the screening of 
unknown mutations of entire human mtDNA were selected for the 
evaluations. The primer set designed by Meierhofer et al., 
containing 48 primer pairs, was originally used for DHPLG 
analysis [7] , whereas the 1 7 primer pairs designed by Bannwarth 
et al. were used in the SN analysis, followed by agarose gel 
detection [11]. These 2 primer sets were designated as the DM 
primer set and SB primer set, respectively, in this study according 
to the names of the first authors in these two publications. 

Materials and Methods 

Ethics Statement 

The parts of this study that involved the use of human blood 
samples from 2 healthy human volunteers and the analysis of 
mtDNA sequences for the blood DNA samples from these 2 
subjects were conducted with approval from the Institutional 
Review Boards of Chang Gung Memorial Hospital in accordance 
with the Declaration of Helsinki for Human Research. The blood 
samples were obtained with written informed consent. 

Cell Lines and Culture of Cells 

143B cells, an osteosarcoma cell line (ATCC CRL 8303), 143B- 
p" cells (clone 87 derived from 143B cells), and a cybrid cell line 
were obtained from Professor Douglas Wallace. The cybrid cells 
line was generated by fusion of the 143B-p'' cells and the cytoplast 
from patients, followed by further selection of a wild-type cybrid 
clone as a control cybrid of the cybrid containing a pathogenic 
mtDNA mutation [23] . 143B cells and the cybrid cells were grown 
in Dulbecco's modified Eagle's medium (DMEM) containing 10% 
fetal bovine serum, 4.5 g/L glucose, and 110 |ig/mL pyruvate. 
143B-p° cells were cultured in the same medium with additional 
supplementation of 50 |J,g/mL uridine [24]. The growth of all cell 
lines was maintained under the conditions of 5% CO2 and 37°C in 
a GO2 incubator. 

DNA Isolation from Whole Blood and Cultured Cells 

For isolation of DNA from whole blood and cultured cells, the 
Centra Puregene Blood Kit and Centra Puregene Cell Kit 
(QIAGEN) were used, respectively, and the procedures described 
in the Cenetra Puregene Handbook (third edition) from QIAGEN 
were followed. The DNA obtained by this way contained both 
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nDNA and mtDNA and so was referred as total DNA in this 
paper. 

Isolation of Mitochondria from Cultured Cells for Isolation 
of mtDNA and the Verification on the Purity of Isolated 
mtDNA 

The procedures of mitochondrial isolation were modified from 
our previous protocols for Western blot analysis [25] with 
additional centrifLigation steps on supematants after the first 
lOOOx^ centrifugation of the homogenate and the DNase I 
treatment step for the mitochondrial fraction because PGR is a 
much more sensitive detection method than Western blot analysis. 
In brief, cells on dishes were washed by phosphate-buffered saline 
(PBS) and then scraped with the mitochondrial isolation solution 
previously described [25]. Cell pellet collected from centrifugation 
at lOOOXjg was homogenized in the mitochondrial isolation 
solution followed by centrifugation at lOOOx^^ for 10 min at 4°C. 
The first supernatant was further centrifuged at 2000 x^, followed 
by centrifugation at 3000 x^ for the second supernatant. The final 
supernatant was centrifuged at 10000 x^ for 15 min at 4°C. The 
pellet or the isolated mitochondrial fraction was washed and 
resuspended in the mitochondrial isolation solution. DNase I 
solution was added to this mitochondrial suspension to a final 
concentration of 10 |.ig/mL and incubated at 37°C for 5 min to 
remove any nDNA attached to the surface of isolated mitochon- 
dria [26]. The mitochondrial suspension was centrifuged and the 
resulting pellet was washed with the mitochondrial isolation 
solution, followed by the mixing of the pellet with the Cell Lysis 
Solution from Centra Puregene Cell Kit. The same DNA isolation 
protocols used for cultured cells were then carried out to extract 
mtDNA from isolated mitochondria, but the Glycogen Solution 
from the kit was added into isopropanol for the DNA precipitation 
step. 

To verify the purity of isolated mtDNA, we compared the 
results of PGR for mtDNA from mitochondria of 1 43B cells, total 
DNA from 143B cells containing both nDNA and mtDNA, and 
DNA from 143B-p° cells containing only nDNA by using nDNA- 
specific primers and mtDNA-spi;cific primers. The mtDNA- 
specific primers that did not amplify DNA from 143B-p cells 
for the amplicon of nt 3374-3894 have been previously described 
[4]. The forward primer 5'-GAAACTTAGATCCAGCTCT- 
CGG-3' and reverse primers 5'-AGCAGCAATTTGTAAGTG- 
TCCG-3' were designed as nDNA-specific primers, which 
amplified a fragment of human S0D2 gene across an exon and 
an intron with the length of 293 bp. 

Automatic Direct Sequencing and Sequence Analysis 

For blood DNA from 2 healthy volunteers (subject A and 
subject B), entire mtDNA was sequenced to obtain information on 
nucleotide positions with SNPs between these 2 subjects. The 
primers from MitoSEQr Resequencing System (ABI) were used to 
generate 46 ampUcons of human mtDNA by PGR according to 
the manufacturer's instruction of the AmpliTag Gold PGR Master 
Mix (ABI). On the other hand, when there were obvious 
heteroduplex signals during DHPLC analysis, the PGR products 
of the DM or SB amplicons from 143B-p" cells, 143B cells, and the 
cybrid cells were sequenced by using the DM or SB primer pairs. 
The PGR products were sent to Genomic Inc. (Taiwan) for 
automatic direct sequencing, which was carried out by the ABI 
3730XL System. Sequencing results were compared with the 
revised Cambridge reference sequence/rGRS (GenBank ID: 
NG_012920) by using the SeqScape software version 2.5 (ABI). 



Cloning of mtDNA Fragments from Blood DNA and 

Subsequent Plasmid Purification 

Based on the sequencing results, we found that there was only 
one SNP at nt 12705 in the amplicon of nt 12566-12895 from 
mtDNA, in which subject A and B had T allele and G allele, 
respectively. This fragment with the size of 330 bp from subject A 
and B was amplified by PGR by using the forward primer 5'- 
GGGTAAGCTTCAAACTAGAGTAC-3' and reverse primer 5'- 
GTAAGGGGAGGATGAAACCGATA and the Taq DNA poly- 
merase (Invitrogen). The PGR products were cloned by using the 
T&A Cloning Vector Kit and the EGOS competent cells 
(Yeastern Biotech, Taiwan) according to the manufacturer's 
instruction. Two sequence-verified colonies with T and C allele 
were amplified. Plasmids purified by the GENO Mini Plasmid Kit 
(Genomics BioScience and Technology, Taiwan) from these two 
clones for allele T and G allele were designated as plasmid A and 
plasmid B, respectively. 

Bioinformatic Analysis to Evaluate the Matching 
Conditions of Primers with NUMTs 

As previously described in our publication [4], to examine how 
primers designed for human mtDNA matched with NUMTs, we 
routinely conducted searches by using the Basic Local Alignment 
Search Tool (BLAST) not only against the nucleotide database in 
National Center for Biotechnology Information (NCBI), but also 
against a local database containing the rCRS and the 228 
accession numbers currently available on GenBank from the 247 
accession numbers that were identified to contain NUMTs 
sequences by Mishmar et al. [13], which were also available on 
the MITOMAP (http://www.mitomap.org/MITOMAP/ 
PseudogeneList). In this study, we used the BioEdit online tool 
(http://www.mbio.ncsu.edu/bioedit/bioedit.html) to create the 
local database. The default settings for the algorithm parameters 
were used for the BLAST searches on either NCBI or BioEdit. 
This BLAST searches were used in this study to examine whether 
any primer pair from DM and SB primer sets had high probability 
to amplify NUMTs. When there were matching between a primer 
and nDNA sequences, the results were categorized into 3 
matching conditions: perfect match, one mismatch, and two or 
more mismatches. When there was one perfect match or only one 
mismatch, the E values from the BLAST searches were all below 
0.3. Moreover, we also examined whether a primer pair matched 
to the same contig on the same chromosome. Subsequendy, the 
nucleotide locations and lengths of the marched regions matched 
by a specific primer pair could be known. 

PCR Reactions for DHPLC and SN/WAVE-HS Analyses 

Optimase DNA polymerase (Transgenomic) and the GeneAmp 
PCR system 9800 (Applied Biosystems) were used for all PGR 
reactions. For the 12566-12895 amplicon, the PGR mixture in 
30 (iL contained the reaction buffer, 1.5 mM MgS04, 0.4 jiM 
primers, 0.2 mM dNTP, 0.75 U polymerase, and 0.03 ng plasmid 
DNA or 6 ng blood DNA. The PCR conditions for this amplicon 
followed the suggestion described in the manual for the SN kit 
except the use of 35 cycles and 64°G as the annealing temperature. 
When using DM primers, PGR was carried out according to the 
conditions described by Meierhofer et Eil., in which one condition 
was applied for all amplicons [7] , except the use of 36 ng DNA per 
30 )iL reaction. For SB primers, the same PCR conditions 
described by Bannwarth et al. [1 1] were used for all amplicons, in 
which 30 ng DNA was added per 30 |4,L reaction. 
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Figure 1. Comparison of the SURVEYOR Mutation Detection Kit 
and the SURVEYOR Plus iVIutation Kit. A mixture of 50% plasmid A 
and 50% plasmid B was used as the template for the PCR. The expected 
size of the PCR product was 330 bp, whereas that of the two DNA 
fragments after SN cleavage was 130 bp and 190 bp. Results for the 
DNA samples without (uncut) and with (cut) SN digestion were 
compared. SN and SN Plus indicate the results obtained using the 
SURVEYOR IVIutation Detection Kit and the SURVEYOR Plus IVIutation Kit, 
respectively. 1/2 SN Plus refers to the use of 50% of the manufacturer's 
recommended amount of Enhancer W2 and Nuclease W when using 
the new kit. The numbers above the peaks of the chromatograms for 
the DNA samples indicate the retention times of the cleaved fragments. 
The numbers above the peaks of the DNA ladder indicate the sizes of 
the DNA markers. The large peak at the retention time after that of the 
300-bp DNA marker indicates the uncleaved DNA. 
doi:1 0.1 371/journal.pone.009281 7.g001 

DHPLC Analysis 

The PCR product amplified fi-om each DNA .sample alone was 
subjected to heating-and-annealing reaction by using the PE2400 
PCR machine (Perkin Elmer) for the generation of homoduplexes 
and heteroduplexes. The samples were first heated to 95 °C for 
4 min, and then cooled to 25°C at the rate of l°C/min. Samples 
were then injected into the WAVE 3500A System equipped with 
the DNASep cartridge (Transgenomic) for repeated analysis of 
each sample at different column temperatures. For the mobile 
phase, buffer A consisted of 0.1 M trietylammonium acetate 
(TEAA), pH 7.0, whereas buffer B contained 25% acetonitrile in 
the TEAA buffer. The gradient profile for each amplicon was 
determined by the Navigator software of the WAVE System. The 
observation of 2 or more peaks on the chromatograms indicated 
the presence of both heteroduplexes and homoduplexes. 

SN/WAVE-HS Analysis 

Heat-annealed PCR product in the amount of 6 |iL was 
subjected to SN digestion by using the Transgenomic SURVEY- 
OR Plus Mutation Detection Kit for Wave and Wave HS Systems 
for most of samples according to the instruction of the 
manufacturer except the use of half amount of Enhancer W2 
and the Nuclease W. For the purpose of comparison, the 
Transgenomic SURVEYOR Mutation Detection Kit, a previous 
version of the kit that has been used in the literature, was also 
tested to perform SN digestion of plasmid DXA. for which 1 of 
Enhancer W and 1 |lL of Nuclease W were mixed with 6 |iL of 
DNA. When making uncut controls, the same DNA samples were 
mixed with all components except the Enhancer and Nuclease, 
which were replaced by water, and subjected to the same reaction 
conditions. DNA samples before (uncut groups) and after SN 
digestion (cut groups) were analyzed simultaneously by the 
ultraviolet (UV) detector (WAVE System) and the fluorescent 



detector (WAVE HS System), which was coupled with a device to 
mix the WAVE HS Staining Solution I (Transgenomic) with DNA 
samples before injection, at the column temperature of 50°C 
under the mode of Double-strand (DS) Multiple Fragments. 
Therefore, the separation of DNA was based on the size of dsDNA 
under a non-denaturing temperature. The chromatographic 
setting was designed according to the instruction in the 
SURVEYOR Plus Mutation Detection Kit when handling 
amplicons from DM primers, which was based on the amplicon 
size, but the minimum and maximum lengths for the analysis were 
set as 100 bp and 1050 bp, respectively, for amplicons from SB 
primers to shorten the run time. A 100-bp DNA ladder (New 
England Biolabs) was run in parallel with samples for each analysis 
under difierent chromatographic conditions. The identification of 
heteroduplexed DNA was determined by the observation of 
additional 2 or more peaks in the chromatograms from a SN- 
digested sample compared with that from the corresponding uncut 
control. 

Results 

Comparison of the Performance of SN/WAVE-HS and 
DHPLC for the Detection of mtDNA Heteroplasmy by 
Using Plasmid DNA and Blood DNA 

So far SN analysis appeared in the literature all utilized the 
SURVEYOR Mutation Detection Kit (old SN kit) from 
Transgenomic, but only the new version of kit containing a new 
Enhancer formula and additional cofactors, SURVEYOR Plus 
Mutation Detection Kit (new SN Plus kit), was available since the 
year of 2010. We therefore compared the results from the 2 
versions of the kits by using the template containing equal amount 
of plasmid A and plasmid B (50% plasmid B). After SN digestion, 
theoretically two fragments in the lengths of 139 bp and 190 bp 
should be generated and detected by SN/WAVE-HS. As shown in 
Figure 1, although we observed that SN digestion performed by 
the new kit, according to the manufacturer's instructions, 
produced higher signals for cleaved products when detected by 
the WAVE HS System, the new kit led to the generation of 
broader peaks, or peaks with shifted retention time, compared 
with those generated by the old kit. This finding was reproducible 
in another 2 independent experiments. We also observed the same 
phenomenon when using a mixture of Control G and Control C 
plasmids provided in the kit (data not shown). These results 
indicated substantial effects of non-specific digestion by SN when 
using the new kit. To resolve such problems without compromis- 
ing the intensities of signals, we used 50% of the manufacturer's 
recommended amounts of Enhancer W2 and Nuclease W 
(Figure 1). 

To compare the sensitivities of DHPLC and SN/WAVE-HS 
analyses for the detection of mtDNA heteroplasmy, we first tested 
the nt 12566-12895 amplicon with the SNP at position 12705 by 
mbdng plasmid A and plasmid B to generate templates containing 
2%-50% of plasmid B. This process excluded any possible 
interference from low-percentage of heteroplasmy undetectable by 
sequencing in blood DNA. The results from DHPLC analysis 
showed that the peaks of the heteroduplexes could not be 
differentiated from the traces of 100% plasmid A or B until at 
least 5"/o of plasmid B remained in the template for the PCR 
(Figure 2A). In contrast, peaks for 2 cleaved products detected by 
SN/WAVE-HS were evident, although small, for the sample 
containing 2% plasmid B, in comparison with the traces of 100% 
plasmid A or B (Figure 2B). This finding indicated that the 
sensitivity of this platform was relatively, although marginally, 
higher than that of DHPLC. As expected, the peaks of the cleaved 
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Figure 2. Evaluation of the detection of low-percentage variants on plasmids by using (A) DHPLC and (B) SN/WAVE-HS. The 

SURVEYOR Plus Mutation Kit was used for SN/WAVE-HS analysis. Templates with different percentages of plasmid B (indicated as plas. B) were 
prepared by mixing plasmid B with plasmid A. The appearance of at least 2 adjacent peaks on the chromatograms from DHPLC analysis indicated the 
presence of heteroduplexes (front peak) and homoduplexes (back peak) in the heat-annealed PGR products. The front heteroduplexes are indicated 
by the symbol of *. Two symbols of T indicate the expected 1 39-bp and 1 90-bp DNA fragments after SN digestion of heat-annealed PGR products. 
The numbers above the peaks of the DNA ladder indicate the sizes of the DNA markers. The small peaks at 3 min were caused by noise from the 
system, which appeared in all chromatograms in the DHPLC analysis; therefore, were disregarded. 
doi:1 0.1 371/journal.pone.009281 7.g002 



products simultaneou.sly detected by the UV detector on the 
WAVE System were considerably less obvious than those detected 
by the flu{)r(;scent detector for the same DNA sample shown in 
Figure 2B because the cleaved products for the 10% plasmid B 
were barely detectable (Figure SI). The traces from the 2% 
plasmid B could not be differentiated from those from pure 
plasmid A or plasmid B when using the old version of SN kit 
(Figure S2), indicating that the relative sensitivity of the old kit was 
slightly lower than that of the new kit. We then applied the 
comparison shown in Figure 2 to samples containing blood DNA 
from subject A and subject B, designated as blood A and blood B, 
respectively, in different percentages. When comparing the 
chromatograms for blood samples containing different percentages 
of blood B widi those for 100% blood A or 100% blood B, the 
peaks of heteroduplexes detected by DHPLC (Figure 3A) and 
cleaved fragments detected by SN/WAVE-HS (Figure 3B) were 
noticeable when the percentage of blood B was > 2%. However, 
100% blood A or lOO'/o blood B (particularly blood B), appeared 
to exhibit a low percentage of heteroplasmy, despite sequencing 
results indicating homoplasmic status. Therefore, comparisons of 
sensitively could be problematic when using blood DNA. 

Evaluation of the Potential of the DM Primer Set for 
Amplifying NUMTs by Bioinformatic Analysis and PCR 
Using the DNA of 143B-p°Cells 

To evaluate whether any primer pair among the DM primer set 
had high potential to amplify NUMTs, we first performed a 
BLAST search for these primers against the NCBI nucleotide 
database and our local database created by BioEdit. We identified 
23 primer pairs among the 48 primer pairs with high possibility of 
matching the sequences of accession numbers from the same 



chromosome under different matching conditions (Table 1). 
Notably, 14 primers pairs had perfect matches with the sequences 
of the accession numbers in one of the databases. Table 1 displays 
the accession numbers of matched nDNA, positions of the 
accession numbers matched by the primers, the sizes of matched 
NUMTs, and information on the regions of mtDNA to be 
amplified. These primer pairs matched with the sequences on 
chromosomes 1, 3, 5, 6, 11, 17, and X, with highest frerjuency for 
chromosomes 1, 5, and 6. The lengths of the highly matched 
NUMTs were the same as those of the targeted amplicons on 
mtDNA. To simpKfy our testing procedures, we then selected 1 7 
primer pairs that had at least one primer that matched perfecdy 
with sequences of at least one accession number of nDNA to 
further confirm the amplification of NUMTs by using the DNA of 
143B-p" cells. The amplicons tested were DM#3, DM#4, 
DM#5, DM#8, DM#12, DM#13, DM#15, DM#16, 
DM#17, DM#18, DM#19, DM#20, DM#21, DM#22, 
DM#24, DM#30, and DM#38. These primer pairs, except 
those for DM#13, DM1#17, and DM#30, exhibited perfect 
matches with sequences of at least one accession number for the 
same chromosome (Table 1). The results from agarose gel 
electrophoresis showed that these selected DM primers could 
produce visible PCR products from the DNA of 143B-p° cells, 
although the yield from 143B-p'' cells was less than that from 143B 
cells, confirming that all the primers selected by our bioinformatic 
approach could all amplify NUMTs (Figure 4A). The PCR 
products from 143B-p" cells with different sizes could be identified 
from the gel for several amplicons, such as DM#8, DM#12, and 
DM#20, although the sizes of major products appeared to be 
similar to those of targeted mtDNA sequences. However, these 
bands, which occurred in different sizes, were not present in the 
lanes for total DNA from 143B cells. Because we observed some 
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Figure 3. Evaluation of the detection of low-percentage variants on Kiuman blood DNA by using (A) DHPLC (B) and SN/WAVE-HS. 

The SURVEYOR Plus Mutation Kit was used for SN/WAVE-HS. Templates with different percentages of the DNA from blood B (designated as bid. B) 
were prepared by mixing the DNA from blood B with the DNA from blood A. The front heteroduplex peaks detected by DHPLC are indicated by *. 
T indicates the expected 39-bp and 190-bp DNA fragments after SN digestion. The numbers above the peaks of the DNA ladder indicated the sizes 
of the DNA markers. 
doi:l 0.1 371/journal.pone.009281 7.g003 



nonspecific products from DM^^H primers, although in low 
quantities, we later excluded this amplicon for furtlier testing as 
obtaining single band should be a common practice before further 
conducting mtDNA detection with the DHPLC or SN/WAVE- 
HS analysis. 

Evaluation of the Potential of the SB Primer Set for 
Amplifying NUMTs by Bioinformatic Analysis and PGR 
with the DNA of 143B-p°Cells 

We further performed the same evaluations on the 17 SB 
primer pairs, which generated long amplicons of approximately 
2 kb. Results showed that 9 primer pairs exhibited high potential 
for matching with NUMTs, although only 2 primer pairs, SB#5 
and SB#13, perfectly matched the same chromosome (Table 2). 
These primer pairs could potentially match sequences on 
chromosomes 1, 6, 8, 17, and X, with highest frequency for 
chromosomes 1, 5, and 6, which was identical to our observation 
in the DM primers. In contrast to our findings in the DM primer 
pairs, the lengths of several NUMTs matched by SB primer pairs 
slightly differed from those of the targeted mtDNA sequences. We 
further selected the primer pairs for SB#2, SB#5, SB#7, SB#8, 
SB#10, SB#12, SB#13, and SB#15 to conduct PGR analyses. 
As shown in Figure 4B, our results confirmed that the selected 
primer pairs could amplify the NUMTs because all of the selected 
SB primers could generate PGR products from the DNA of 143B- 
p" cells. The sizes of these PGR products were very close to those 
from 143B cells, although the amounts of the PGR products were 
relatively lower for amplicons SB^5, SB^S, and SB^15. 

Demonstration on the Purity of mtDNA Isolated from the 
Mitochondria of Cultured Cells 

To clearly differentiate true NUMT effects from real mtDNA 
heteroplasmy or possible PGR artifacts, it was essential to compare 
the results from total DNA isolated from whole cells with those 



from pure mtDNA isolated from mitochondrial fraction for 
suspected primer pairs. To achieve this objective, it was necessary 
to obtain mtDNA without nDNA contamination, and to verify the 
purity of the isolated mtDNA. As shown in Figure 5, the primer 
pair for the 3374—3389 amplicon of mtDNA was mtDNA-specific 
because it could amplify targeted mtDNA sequences from the total 
DNA of 143B cells or cybrid cells, or the mtDNA of isolated 
mitochondria, but not the DNA of 143B-p" cells. When we used 
the primer pair for the nDNA-encoded S0D2 gene, the DNA of 
143B-p'' cells, 143B cells, or cybrid cells could be amplified, 
whereas the DNA from isolated mitochondria could not. These 
examinations proved that the mtDNA purified using our 
procedures did not have nDNA contamination. 

Assessment of the Effects of NUMTs from the DM Primer 
Set on mtDNA Detection during DHPLC and SN/WAVE- 
HS Analysis 

We evaluated the DM primer pairs displayed in Figure 4A, 
except DM#8, by comparing the results for total DNA and those 
for isolated mtDNA of the same cells analyzed in the same rum. 
Each DNA sample was subjected to heat-annealing reactions. 
Among the 16 aiuplicons tested, 8 amplicons displayed small, but 
clear, heteroduplexes on the chromatograms of DHPLC analysis 
at certain temperatures for 143B cells, whereas 4 amphcons 
displayed ambiguous shoulder peaks adjacent to the major 
homoduplexes. However, we observed identical patterns for the 
mtDNA of 143B cells (chromatographic data not shown). Similar 
results could be found when using 143B cybrid cells, although the 
optimal analyzing temperatures and the amplicons with clear 
heteroduplex peaks differed for a few ampUcons because of 
different mtDNA background (chromatographic data not shown). 
However, sequencing results did not reveal any obvious hetero- 
plasmic sites for these amplicons with heteroduplex signal detected 
during DHPLC analysis (data not shown). Therefore, we could not 
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Figure 4. PCR products amplified from the DNA of 143B and 143B-p° ceils by using the selected primer pairs. (A) The 17 primer pairs 
selected from the DM primer set. (B) The 8 primer pairs selected from the SB primer set. A 100-bp DNA ladder and a 1-kb DNA ladder were used as the 
DNA sizing markers for the testing of the DM primer pairs and SB primer pairs, respectively. N, no-template control; R, DNA of MBB-p" cells as the 
template; B, DNA of 143B cells as the template. 
doi:1 0.1 371/journal.pone.009281 7.g004 



confirm the interference of NUMTs on the interpretation of the 
results of DHPLC analysis for these ampUcons and primer pairs. 
Table 3 lists the results from these evaluations. 

The same evaluation was then applied for SN/^\'AVE-HS 
analysis. It is important to emphasize that the SN/^VAVE-HS 
analysis not only can detect heteroduplexes and predict approx- 
imate cleavage site by comparing the patterns prior to (uncut) and 
after (cut) SN digestion, but also can detect PCR products of 
different sizes. Our results for the total DNA of 143B cells showed 
that among the 16 DM amplicons analyzed, 5 amplicons 
displayed clear peaks other than the major uncut peaks in uncut 
groups and also obvious heteroduplexes, indicated by the changes 
in the patterns of peaks before and after SN digestion. However, 
except for DM^38, the same peaks or patterns of changes were 
also found when analyzing the mtDNA of 143B cells (chromato- 
graphic data not shown). We also analyzed these 5 amplicons from 
the cybrid cells. Except for DM#38, the chromatograms for the 
total DNA and the mtDNA of cybrids cells were also identical for 
each amplicon, but the patterns of several amplicons were different 
from those of 143B cells (chromatographic data not shown). 
Table 3 lists the summary of these findings and Figure 6A displays 
the chromatograms for DM^38, which was a reprehensive result 
from 2 repeated experiments. We observed one peak <100 bp in 
the uncut groups in the total DNA, but not the mtDNA, for 
amplicon DM#38 of the 143B and cybrid cells. The peak height 
of another large peak of 1 00-200 bp for the uncut groups of the 
total DNA was decreased to near the background level of the 
uncut groups of the mtDNA from both types of cells. Moreover, 



after SN digestion, the shift of the 100-200 bp peak to a smaller 
peak for the total DNA of the 143B and cybrid cells did not occur 
for the mtDNA. These results demonstrated the interference from 
NUMTs on SN/WAVE-HS analysis of DM#38 in both the 143B 
and cybrid cells. In contrast, as shown in Figure 6B, the differences 
between the total DNA and the mtDNA in DHPLC analysis of 
143B cells or the cybrid cells were absent, although the 
chromatograms for the total DNA and the mtDNA both displayed 
ambiguous shoulder peaks adjacent to the homoduplex peaks. 

Assessment of the Effects of NUMTs from the SB Primer 
Set on mtDNA Detection during SN/WAVE-HS Analysis 

The same strategy was further applied to evaluate selected SB 
primer pairs during SN/WAVE-HS analysis. It was not tested for 
DHPLC analysis because SB amplicons were too long to be 
analyzed by DHPLC. After SN digc-stion, all tested amplicons 
produced several peaks in low abundance in the total DNA of the 
143B or cybrid cells, but the patterns were highly similar to those 
of the mtDNA (data not shown) except SB#12. A representative 
data from 2 repeated experiments were displayed in Figure 7. As 
shown in the upper panel of Figure 7, there was a dominant peak 
of approximately 800 bp in the uncut groups for the total DNA 
and the mtDNA from 143B cells, but 2 peaks at the similar 
retention time showed up in the cut group of the total DNA, but 
not the mtDNA. On the contrary, we observed a peak close to 
1000 bp in the chromatogram for the mtDNA, but not the total 
DNA, of 143B cells. These findings indicated the contributions of 
NUMTs to such differences. However, the results displayed in the 
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Figure 5. Confirmation of tlie purity of isolated mtDNA from 143B cells and the cybrid cells. The purity of isolated mtDNA was 
determined by comparing the PCR results obtained using mtDNA-speclfic and nDNA-specific primers for DNA from isolated mitochondria, 143B-p° 
cells, 143B cells, and the cybrid cells. The primer pair for the amplicon 337439-bp and 190-bp 3894 of human mtDNA (HIVI3374-3894) was used as the 
mtDNA-speclfic primers, which generated PCR products with an expected size of 521 bp. One primer pair for the S0D2 gene was used as the nDNA- 
specific primers, which generated PCR products with an expected size of 290 bp. N, no-template control; p°, DNA of 143B-p° cells as the template; T, 
total DNA from whole cells as the template; m, mtDNA extracted from isolated mitochondria as the template. 
doi:1 0.1 371 /journal.pone.009281 7.g005 



lower panel of Figure 7 indicated that the patterns of peaks in the 
cut groups did not exhibit obvious differences between the total 
DNA and the mtDNA of the cybrid cells. 

Sequencing Results for the DM#38 and 5B#12 
Amplicons of 143B-p°Cells, MSB Cells, and the Cybrid 
Cells 

To examine whether sequencing data of different cells could be 
used to explain the findings for the DM^38 and SB7^12 
amplicons during SN/WAVE-HS analysis, PCR products for 
these 2 amplicons from 143B-p" cells, 143B cells, and the cybrid 



cells were sequenced. The differences in the sequences among 
these 3 ceU lines for the nt 12806-13311 (DM#38) and nt 4995- 
7048 (SB#12) regions were Ksted in Table SI and Table S2, 
respectively. Our sequencing results showed that the mtDNA of 
the 143B cells and cybrid cells differed only at position 13135 for 
DM#38, and that mtDNA of these two cell lines had 14 
nucleotides with the sequences different from those of 143B-p'' 
cells (Table SI). However, although the sequences of the DNA of 
143B-p" cells also differed from mtDNA sequences of either of the 
143B or cybrids at 26 nucleotide positions for SB#12, 6 of 7 
nucleotides positions with different mtDNA sequences between 
143B cells and the cybrid cells had the same sequences between 



Table 3. Summary on the evaluation of DM primers during DHPLC and SN/WAVE-HS analyses. 





DHPLC 




SN/WAVE-HS 








143B 


Cybrid 


143B 




Cybrid 




DM amplicon 






Uncut 


Cut vs. Uncut 


Uncut 


Cut vs. Uncut 


#3 


Y 


Y 


N 


N 


ND 


ND 


#4 


Y 


Y 


Y 


Y 


Y 


Y 


#5 


Y 


Y 


N 


N 


ND 


ND 


#12 


N 


ND 


N 


N 


ND 


ND 


#13 


7 


Y 


N 


N 


ND 


ND 


#15 


Y 


Y 


Y 


Y 


Y 


Y 


#16 


Y 


7 


Y 


Y 


Y 


Y 


#17 


Y 


7 


N 


N 


ND 


ND 


#18 


N 


ND 


N 


N 


ND 


ND 


#19 


? 


7 


N 


N 


ND 


ND 


#20 


Y 


Y 


N 


N 


ND 


ND 


#21 


N 


ND 


N 


N 


ND 


ND 


#22 


Y 


Y 


Y 


Y 


Y 


Y 


#24 


? 


Y 


N 


N 


ND 


ND 


#30 


N 


ND 


N 


N 


ND 


ND 


#38 


? 


7 


Y 


Y 


Y 


Y 



For DHPLC analysis, Y indicates the presence of small but clear heteroduplex peaks, whereas N indicates the absence of such peaks. The question mark indicates the 
presence of ambiguous shoulder peaks. For SN/WAVE-HS analysis, Y in the uncut groups indicates the presence of cleaved DNA fragments of sizes different from those 
of the undigested peak, whereas Y in the cut vs. uncut column indicates changes in patterns of peaks from cleaved fragments after SN digestion. ND indicates that the 
experiments were not performed for cybrid cells because there was no positive signal or change in the chromatograms of 143B cells. 
doi:1 0.1 371 /journal.pone.009281 7.t003 
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Figure 6. Results from SN/WAVE-HS and DHPLC analysis of the DM#38 amplicon by using total DNA and mtDNA from cells. (A) 

Results from SN/WAVE-HS analysis. (B) Results from DHPLC analysis. The upper and lower panels of Figure 6A and Figure 6B display the results for 
143B cells and the cybrid cells, respectively. In Figure 6A, the results for the DNA samples without (uncut) and with (cut) SN digestion were also 
compared. T indicates the peaks that appeared in the chromatograms of uncut DNA but were not observed in the chromatograms of the 
corresponding cut DNA for all DNA samples. V indicates the peaks >100 bp only present in the chromatograms of the uncut DNA of both total DNA 
and mtDNA (m), but not the cut DNA. Such peaks for the mtDNA were considerably smaller than those of total DNA. I indicates the peaks <100 bp 
that were present in the chromatograms of the uncut total DNA but absent in those of the uncut mtDNA. The patterns of peaks in the 
chromatograms from uncut DNA and cut DNA displayed marked differences. The chromatograms of the cut DNA samples displayed numerous small 
peaks. One obvious peak (indicated by *) appeared for the cut DNA of whole cells but was absent for mtDNA. In Figure 6B, the small but identifiable 
heteroduplex peaks detected by DHPLC are indicated by #. 
doi:1 0.1 371/journal.pone.009281 7.g006 



the mtDNA of the cybrid cells and the DNA of 14,3B-p cells, but 
not between the mtDNA of 143B ceUs and the DNA of 143B-p" 
cells for this long amplicon (Table S2). Our sequencing results 
might explain the identical NUMT effects observed in 143B cells 
and the cybrids cells for DM#38, because of one nucleotide 
difference in the mtDNA background between these 2 cell lines for 
this region. The clear NUMT effect in 143B cells, but not the 
cybrid cells, for SB#12 might be caused by closer sequence 
similarity between nDNA of 143B-p° cells and mtDNA of the 
cybrid cells than between 143B-p° cells and 143B cells. However, 
the chromatographic data of sequencing indicated that the 
mtDNA sequences in these 2 regions for 143B cells and the 
cybrid cells displayed no sign of heteroplasmy (data not shown). It 
was therefore impossible to identify the presence of interference 
from NUMTs based on the sequencing results for the DM#38 
and SB#12 ampHcons. 



Discussion 

Our current study is the first to apply the SN/WAVE-HS 
platform for the detection of mtDNA and to compare its 
performance with that of DHPLC. Our results suggest that the 
relative sensitivity of the SN/WAVE-HS platform for the 
detection of mtDNA heteroplasmy is marginally higher than that 
of DHPLC for one SNP, although both techniques determine 
heteroplasmy qualitatively. Furthermore, we have provided a 
practical approach for the screening and evaluation of the 
potential interference of NUMTs on mtDNA detection during 
DHPLC and SN/WAVE-HS analysis from primers by using 
published primer sets as an example. Most importantly, this study 
is the first to verify NUMTs effects on the determination of 
mtDNA heteroplasmy by obtaining pure mtDNA from isolated 
mitochondria without PCR-detectable nDNA contamination. 

Both DHPLC and SN/WAVE-HS platforms are qualitative 
techniques and thus was not suitable to be evaluated for the exact 
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Figure 7. Results from SN/WAVE-HS analysis of the SB#12 
amplicon by using total DNA and mtDNA from cells. DNA 

samples from isolated mitochondria (m) and whole cells of 143B cells 
and the cybrid were used. The upper and lower panels display the 
results for 143B cells and the cybrid cells, respectively. The results for 
the DNA samples without (uncut) and with (cut) SN digestion were also 
compared. T indicates the peaks that appeared in the chromatograms 
of uncut DNA but were absent in the chromatograms of the 
corresponding cut DNA for all DNA samples. Two obvious peaks 
(indicated by *) showed up after SN digestion for the total DNA of 143B 
cells, but could not be apparently visualized for the mtDNA of 143B 
cells (upper panel). I indicates the peak close to 1000 bp in the cut 
DNA of the mtDNA, but not in the cut DNA of the total DNA from 143B 
cells (upper panel). These two differences in the patterns of peaks 
between total DNA and mtDNA found in 143B cells were absent in the 
cybrid cells (lower panel). 
doi:l 0.1 371/journal.pone.009281 7.g007 

limit of detection. By conducting analysis on a serial of samples 
with different percentages of mutant, Janne et al. reported that the 
sensitivity of SN/WAVE-HS method for the EFRF gene reached 
1 % of mutant [10], which was similar to 2% observed in this study 
for plasmids with an SNP. We identified that DHPLC was suitable 
to detect mutant at 5% or greater in this study, which was 
consistent with that reported in our previous study [4] and in other 
studies [8,27] for mtDNA. Because testing on the optimal 
temperatures required for DHPLC analysis is not conducted for 



SN/WAVE-HS analysis, and the WAVE HS System is superior to 
gel electrophoresis for the identification of cleaved DNA fragments 
with their approximate sizes, theoretically, the SN/WAVE-HS 
platform should provide a considerably more efficient and 
convenient method than DHPLC for the screening of mtDNA 
mutations and heteroplasmy. However, we only tested the template 
with one mismatch in the ampUcon. When there are more than one 
mismatch in one amplicon, the signals indicating the presence of 
heteroduplexcs and sensitivity should be enhanced during DHPLC 
analysis, but might be reduced in SN/WAVE-HS analysis due to 
the formation of many small cieaved fragments, especially that the 
size of resulting cleaved DNA fragments were recommended to be 
> 100 bp for the detection in the WAVE HS System. On die other 
hand, our results suggest that several factors, including non-specific 
activities of SN, the presence of co-amplified NUMTs, the presence 
of non-specific PCR products of different siz(-s prior to SN digestion, 
potential PCR artifacts, and true low-abundance of mtDNA 
heteroplasmy that commonly occurs in samples, may cause highly 
ambiguous signals during SN/WAVE-HS analysis. Such problems 
may be amplified when working on large amplicons. Therefore, 
although we did not test amplicons with multiple SNPs from 
mtDNA of human blood for the SN/WAVE-HS analysis in this 
study, we suspected that the usefulness of the SN/WAVE-HS 
platform in the screening or identification of mutations for mtDNA 
may not be as good as for nDNA. 

Although our methods consisting of bioinformatic analysis and 
PCR verification using the DNA from HSB-p" cells indicated the 
amplification of NUMTs by several primer pairs, only 2 of 24 
primer pairs tested experientxd interference by NUMTs during 
SN/WAVE-HS analysis. None of these primers were associated 
with NUMT effects during DHPLC analysis after comparing the 
results obtained using pure mtDNA and total DNA for either 143B 
cells or the cybrid cells. Therefore, low-level heteroduplexes, 
which could be regarded as low-level heteroplasmy of mtDNA, 
detected in some amplicons from 143B cells could be true 
heteroplasmy at any positions in these amplicons. They could be 
also heteroduplexes and chimerical molecules generated during 
PCR [28,29]. One critical implication of our results was that the 
use of mtDNA deriving from isolated mitochondria without 
nDNA contamination was essential to prove the interference from 
NUMTs for a primer pair. However, so far only the study of 
Calvignac et al. has investigated the feasibifity of extracting 
mtDNA from mitochondrial fraction from the superfamily 
Aselloidea on preventing contamination of co-amplified NUMTs 
during PCR, but they simply conducted one low-speed centrifu- 
gation followed by one high-speed centrifugation before the DNA 
extraction of the final pellet without showing any data to 
demonstrate the purity of isolated mtDNA [22]. In this study, 
we found that one low-speed centrifugation was not sufficient to 
fully centrifuge down all nuclear fractions and cell debris, and 
treatment of the final mitodiondrial pellet by DNase I was also 
critical to eliminate any nDNA attached to the surface of 
unbroken mitochondria. Moreover, if proper primers and control 
DNA samples were not used to confirm the purity of isolated 
mtDNA as we have performed, the conclusion from the results 
might be misleading. However, our results did not suggest that 
pure mtDNA from real samples should be obtained during routine 
analysis of mtDNA mutations when using high-sensitivity 
techniques because it is not practical for most studies and is 
generally not applicable in clinical studies. Instead, our study 
developed a strategy for the evaluation of NUMT effects during 
primer design by obtaining an "mtDNA control" from unlimited 
cultured cell lines for confirming NUMT effects before using these 
primers to analyze real specimen. Our strategy is useful because p° 
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cells are not available to most investigators, and primer pairs that 
can amplify the DNA from p" cells, which is often hard to be 
avoided for some mtDNA regions, such as COI and COII genes 
with long sequences of high similarity to those of NUMTs, are 
potentially still usable through the evaluation by using our strateg)-. 

Investigators have created databases for NUMTs by mate hing 
entire mtDNA sequences with different human genome databases 
by applying different criteria [13,17], and by sequencing cloned 
PGR product from p" cells [30] or isolated nDNA from sperm cells 
[31] with the use of2 several primer pairs covering all regions of 
mtDNA. Although those databases may be useful for the screening 
of primers matchable to NUMTs, the real probability of co- 
amplification of NUMTs with mtDNA during the PGR by a 
primer pair and the subsequent effects of NUMTs on the detection 
of mtDNA variants might not be predictable by using these 
databases. The reasons are that primer pairs other than the 
perfecdy matched primer pairs are still able to amplify NUMTs, 
and that the abundances of co-amplified NUMTs are dependent 
on the copy numbers of NUMTs amplified by each primer pair. 
For example, in our study, although two primers for the amplicon 
DM#38, which was associated with the NUMT effects, had 
perfect matches with the same chromosomes, only the reverse 
primer for the amplicon SB7^12, also associated with the NUMT 
effects exhibited perfect matches with sequences of multiple 
accession numbers. However, 16 of 20 bp of the forward primer 
for SB#12 in fact matched with at least 3 sites in 2 chromosomes 
without any gap as the reverse primer (data not shown). Another 
problem is that the NUMT databas(;s producxxl by the sequencing 
of cloned PGR products amplified from the DNA of p" cells or 
isolated nDNA by using a set of specific primer pairs might not 
predict the probability of the co-amplification of NUMTs by a 
different primer pair. Moreover, the absence of mtDNA in isolated 
nDNA for such purposes were not demonstrated, such as in the 
study of Ramos et al., although the authors claimed the purity of 
the nDNA used [31]. Therefore, the "m dim primer validation" 
approach employed by Ramos et al. [17], which used a database 
generated from bioinformatics analysis, to make conclusion on the 
risk of the co-amplification of NUMTs with mtDNA, and to 
determine whether the reported mtDNA mutations were authentic 
for the primer pairs used in cancer studies by other studies might 
mistakenly exclude many true mtDNA mutations reported 
previously. 

When heteroplasmic variants of mtDNA were found to match 
with nucleotide sequences in NUMT databases, many investiga- 
tors often inferred that such findings should be caused by co- 
amplification of NUMTs. However, without appropriate control 
experiments, particularly when conducting low-sensitivity se- 
quencing, the conclusion could be considered arbitrary. For 
example, in the study of Parr et al. [30], when PGR products on 
the agarose gel could be generated from the DNA of both p" cells 
and clinical samples with the same primer pairs, a proportion of 
the heteroplasmic signals from sequencing were determined to be 
false results because some identical nucleotide sequences existed in 
the additional variant and the NUMTs. However, no control 
DNA, such as the pure mtDNA employed in this study, was tested 
to verify such conclusion. If nucleotide signals from NUMTs could 
be so easily detected by the low-sensitivity sequencing technique, 
as long as the same primer pairs could amplify the DNA of p" cells, 
the interference of NUMTs to the detection of mtDNA 
heteroplasmy during DHPLG or SN/WAVE-HS analyses would 
have been tremendous. However, our results showed the 
probability of such interference was low, despite several primers 
pairs perfecdy matching the nDNA sequences on the same 
chromosomes. 



NUMTs may affect the detection of mtDNA during DHPLG 
and SN/WAVE-HS analyses in manners additional to the factor 
of NUMT copy number. Heteroduplex formation between 
NUMTs and mtDNA, or between different NUMTs, can 
potentially cause low-level heteroduplex signals during DHPLG 
analysis and the generation of cleaved DNA fragments after SN 
digestion. In contrast, signals of heteroduplexes from true mtDNA 
heteroplasmy might be suppressed because of a dilution effect 
caused by NUMTs containing the same nucleotide sequences as 
those of the dominant allck's. Theoretically, the greater the extent 
of mismatch in one amplicon, the more effective the resolution of 
the heteroduplex signals during DHPLG analysis, whereas it may 
make cleaved DNA fragments too small to be detectable during 
SN/WAVE-HS analysis. Therefore, although we found that low- 
level heteroduplexes detected by DHPLG for some primer pairs 
tested in this study were not related to NUMTs, the likehood of 
generating heteroduplex signals caused by NUMTs during 
DHPLG analysis when using other untested primer pairs remains 
high and requires further evaluation. On the other hand, a unique 
interference from NUMTs that can occur during SN/WAVE-HS 
analysis, but not during DHPLG analysis, is the detection of 
putative NUMT products with sizes differing from those of the 
expected mtDNA amplicons in samples prior to SN digestion. 
However, non-specific PGR products may also cause these results 
because the WAVE HS System can detect DNA fragments in 
different sizes with much greater sensitivity and clarity than 
agarose gel electrophoresis can. Accordingly, the use of isolated 
mtDNA is important to clarify this issue. As shown by the results of 
Figure 6 and Figure 7, peaks present in the chromatograms of both 
total DNA and mtDNA should be caused by non-specific PGR 
products. When the patterns of certain peaks were different 
between the chromatograms of total DNA and mtDNA, NUMTs 
are likely to play the major roles in such differences. Therefore, the 
absence of peaks in the uncut and cut samples of mtDNA in 
comparison with the pattern of total DNA, for DM7^i^:38 indicates 
enhancement of non-specific signals caused by NUMTs (Figure 6). 
On the other hand, our results suggest that samples with different 
mtDNA backgrounds are subject to different types of NUMT 
effects. Gonsequently, we observed potential interference from 
NUMTs in both 143B and cybrid cells for the DM#38 amphcon 
(Figure 6), but only in 143B cells for the SB#12 amphcon 
(Figure 7). It was possibly because of more differences in the 
nucleotide sequences between mtDNA of 143B and mtDNA of the 
cybrid cells, or higher sequence similarity between nDNA of 143B- 
p cells and mtDNA of the cybrid cells than between 143B-p°ceUs 
and mtDNA of 143B ceUs, for the SB#12 amplicon. 

Overall, our studies have produced several crucial novel 
findings. First, the sensitivity of SN/WAVE-HS analysis for the 
detection of mtDNA heteroplasmy is marginally greater than that 
of DHPLG analysis for amplicons with one SNP, but several 
potential factors may affect decision-making on the heteroplasmy 
or homoplasmy status of mtDNA during SN/WAVE-HS analysis. 
Second, we results suggest that the use of a bioinformatics 
approach or the use of the DNA from p" cells for the evaluation on 
the potential amplification of NUMTs by a primer pair is helpful 
during primer design to exclude some primers with high risk to 
amplify NUMTs, but is not capable of fuUy predicting potential 
interference by NUMTs on the judgment of mtDNA heteroplasmy 
during the detection by a high-sensitivity and PCR-based 
technique. Third, mtDNA isolated from the mitochondria of cell 
lines with confirmed absence of nDNA contamination by PGR is 
required during primer design to verify NUMT effects from a 
highly suspected primer pair for subsequent analysis of mtDNA in 
the total DNA from real samples. Finally, when NUMTs interfere 
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with the judgment of mtDNA hrteroplasmy, the extent of such 
effects may be dependent on the mtDNA backgrounds of samples 
and the detection methods used. 

Supporting Information 

Figure SI Simultaneous detection of the same DNA 
samples used in Figure 2B by using the UV detector 
during SN/WAVE-HS analysis. During the detection of the 
DNA by the HS detector on the WAVE System, the data acquired 
by the UV detector of the WAVE System were recorded 
simultaneously. T indicates the expected DNA fragments in the 
sizes of 139 bp and 190 bp after SN digestion of heat-annealed 
PGR products. 
(TIFF) 

Figure S2 SN/WAVE-HS analysis of the same PGR 
products used in Figure 2 by using the old SN kit. The 

same heat-annealed PGR products were analyzed by SN/WAVE- 
HS by using the SURVEYOR Mutation Detection Kit. 
(TIF) 
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